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Abstract 


Standard  perturbation  results  are  applied  to  the  case  of  laser 
induced  dynamic  Stark  shifts  of  vibration-rotation  terms  of  diatomic 
molecules.  For  fields  at  frequencies  near  an  allowed  vibration- 
rotation  transition  there  are  just  four  terms  coupled  by  allowed 
transitions  which  cause  significant  shifts.  After  discussion  of  the 
method  of  calculating  specific  matrix  elements  from  rotationless 
Einstein  A coefficients,  examples  are  given  of  Stark  shift  predictions 
for  terms  of  HF,  DF  and  CO  involved  in  strong  laser  lines  of  these 
molecules. 
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PREFACE 

This  report  contains  essentially  all  of  the  results  obtained 
under  this  contract.  We  were  able  to  establish  the  specific  form  of 
the  laser-induced  Stark  shifts  in  diatomic  molecular  terms,  to  locate 
'*  suitably  reduced  intensity  data  (Einstein  coefficients)  for  the 

originally  proposed  molecules  (to  obtain  the  necessary  dipole  matrix 
elements),  and  to  make  specific  calculations  of  the  frequency  depen- 
dence and  magnitude  of  the  effect  in  these  molecules  (HF,  DF,  and 
CO). 

In  addition  to  the  principal  investigator  the  following 
personnel  were  associated  with  the  project: 

W.  M.  Benesch,  Institute  for  Physical  Science  and  Technology, 
University  of  Maryland 

T.  D.  Wilkerson,  Institute  for  Physical  Science  and  Technology, 
University  of  Maryland 

R.  S.  Klingbeil,  Systems  Planning  and  Analysis,  Inc. 

A trimmed  version  of  the  present  report  is  in  preparation  for 
the  Journal  of  Quantitative  Spectroscopy  and  Radiant  Transfer.  An 
initial  report  of  the  results  was  presented  as  paper  WlA.l  at  the 
International  Conference  on  Multiphoton  Processes,  University  of 
Rochester,  June  1977. 


I.  Introduction 


One  of  the  consequences  of  the  availability  of  tunable  lasers  has 
been  the  revival  of  interest  in  the  subject  of  resonance  fluorescence.  With 
this  have  come  numerous  theoretical  calculations,  notably  Hollow  (1969), 
showing  that  the  emission  and  absorption  spectra  lines  are  split  by  the 
presence  of  an  intense  monochromatic  radiation  field.  The  splitting  on 
resonance  is  found  to  be 


+ n = + (y^j'^)/(hc)  (cm 


(1) 


with  E the  incident  laser  field  amplitude  and  the  dipole  matrix 

element  of  the  transition.  This  splitting  can  be  identified  with  the  Rabl 
frequency  of  population  inversion  (Rabi,  1937)  and  the  "dynamical  Stark  shift" 
of  Autler  and  Townes  (1955).  Very  beautiful  confirmation  of  the  calculations 
has  been  demonstratea  with  Na  atomic  florescence  by  Schuda,  Stroud  and  Hercher 
(1974),  Wu,  Grove  and  Ezekiel  (Wu,  (1975),  Grove,  (1977))  and  Walther's 
group  (Hartig,  1976). 

In  earlier  work  Russian  investigators  had  anticipated  the  importance  of 

laser  induced  Stark  shifts, even  for  non-resonance  cases,  due  to  the  high  available 

fluxes.  The  subject  is  nicely  reviewed  by  Eonch-Bruevlch  and  Khodovoi  (1967). 

They  continued  their  work  and  demonstrated  large  (.1  to  2.5  cm  ^)  shifts  of 

the  potassium  atomic  D fluoresence  line  for  Ruby  illumination  of  30  to  300 
2 

MW/cm  intensity  (Bonch-Bruevich,  1969).  There  is  even  a quantum  electronics 
textbook  description  of  the  effect  with  order-of-magnitude  estimates  in  Pantell 
and  Puthoff  (1969). 

More  recently  the  powerful  output  of  a COj  TEA  laser  has  been  used  to 

produce  large  non-resonant  Stark  shifts  in  the  excited  states  of  atomic  He 

o o 

(Prosnitz,  1976).  Shifts  of  2 to  4 A at  4713  A were  observed  for  incident 
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intensities  having  amplitudes  of  130  to  215  kV/cm. 

Much  less  work  has  been  done  for  case  of  laser  induced  Stark  shifts  in 

molecules,  despite  the  fact  that  the  most  powerful  lasers  themselves  run  on 

molecular  transitions.  Braun  and  Petelln  (Braun,  1974)  calculated  Stark 

shifts  for  CO2  Irradiated  diatomic  molecule  rotation  and  vibration-rotation 

states.  They  estimated  splitting  of  the  lowest  vibrational  levels  of  HCl 

of  .1  cm  ^ and  of  CO  of  .2  cm  ^ for  a CO2  laser  running  at  9.9  microns  with 

6 2 

a field  amplitude  of  10  V/cm  (1.3  GW/cm  ).  These  shifts  are  quite  signifi- 
cant, considering  the  narrowness  of  the  transitions  themselves,  but  are  not 
as  large  as  will  occur  at  resonance,  a case  not  considered  by  Braun  and 
Petelin.  Fedorov  and  Makarov  (Fedorov,  1975)  have  published  a short  note  on 
the  resonance  case  in  diatomic  molecules  but  the  form  of  their  results  is 
not  easily  used. 

Aldridge  et  al.  (1976)  have  attempted  to  estimate  dynamical  Stark  effects 
in  their  studies  of  laser  isotope  separation.  They  used  the  results  of  Braun 
and  Petelin  for  the  non-resonant  case  and  attempted  to  use  the  ideas  of  power 
broadening  to  estimate  the  resonant  splitting  under  the  mistaken  assumption 
that  the  appropriate  Rabi  frequency  is  given  by  the  static  dipole  moment. 

This  leads  to  a large  overestimate  for  HCl. 

There  has  also  been  an  extensive  effort  to  calculate  the  effect  of  laser 
fields  on  the  emission  and  absorption  of  a gas  such  as  a laser  medium.  This 
work  is  reviewed  by  Beterov  and  Chebataev  (Beterov,  1975)  who  point  out  the 
parts  of  the  calculated  lineshapes  due  to  dynamical  Stark  effect.  Particularly 
notable  is  the  paper  of  Skribanowitz,  Kelly  and  Feld  (1972)  who  show  calcu- 
lated examples  of  the  M (magnetic  quantum  number)  dependence  of  the  splitting 
in  molecules.  Typically  these  calculations  are  for  fields  such  that  the 
splittings  are  not  much  larger  than  the  Doppler  llnewldths  and  so  are  aimed 
at  laser  Induced  line-narrowing  spectra  (hole  burning) . 

I 
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The  purpose  of  this  report  Is  to  give  the  form  of  the  dominant  dynamical 
Stark  shift  on  diatomic  molecule  vibration-rotation  terms  in  conventional 
molecular  notation  and  present  some  explicit  calculated  examples  of  these 
shifts  as  a function  of  incident  laser  frequency  as  the  frequency  goes  through 
resonance.  Contact  is  made  with  the  2x2,  degenerate  perturbation 
results  used  in  the  atomic  resonance  case  and  they  are  applied  to  molecules, 
attempt  to  treat  very  carefully  the  calculation  of  the  appropriate  dipole 
transition  matrix  elements  from  tabulations  of  the  Einstein  A coefficients 
of  the  molecular  transitions.  The  examples  chosen  correspond  to  high  power 
lines  of  HF,  DF  and  CO  lasers  where  large  Stark  shifts  should  occur  but 
as  yet  have  not  been  reported. 


II.  Theoretical  Background 

The  conventional  form  of  the  off-resonance  level  shift  of  the  nth 
level  for  an  applied  electric  field  e = cos(kx-wt)  can  be  written 
(Bonch-Bruevich,  1967), 


AE^(u)) 


4 4 ] E -E  -hoj 
m n m 


+ 


E -E  +hu 
n m 


2 


} 


(2) 


with  E the  energy  of  the  mth  level,  p the  dipole  transition  matrix 
m mn 

element  between  levels  m and  n.  Because  of  the  tradition  in  molecular 
spectroscopy  of  giving  the  spectral  constants  in  "wave  number"  units  (cm 
it  is  convenient  to  rewrite  this  as  the  shift  of  the  n^  term, 

6v^  = AE^/(hc).  It  is  also  useful  to  define  transition  Rabl  frequencies 
in  cm  ^ units. 


mn 


•Eg  I /(he) 


(3) 


We 
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Eq.  (2)  can  be  rewritten  as 


6v 

n 


(4) 


with  the  energy  of  the  nth  term  and  the  frequency  of  the  applied 

radiation  field  v,  both  in  cm  Combining  the  two  terms  in  the  curly 
brackets  in  Eq.  (4)  gives 


^ 1 r ^mn  ^'^n  '^m^ 

;f-2 

m (v^-v^)  -V 


(5) 


The  functional  form  in  Eq.  (5)  shows  directly  the  specific  character  of  the 
signs  of  each  contribution  to  the  n^^  term  shift.  This  sign  depends  on 
whether  the  driving  field  frequency  is  smaller  or  larger  than  the  transition 
frequency  magnitude  whether  is  below  or  above  . 

For  terms  spaced  widely  compared  to  the  Rabi  frequencies  , and  v 

not  closer  than  5 to  the  transition  frequency 

accurate  expression  of  the  term  shifts.  We  will  see  below  that  the  allowed 

transitions  to  a given  term  have  frequencies  differing  by  20  cm  ^ or  more 

while  the  Rabi  frequencies  are  less  than  1 cm  ^ for  Intensities  less  than 
8 2 

10  Watts/cm  so  that  near  exact  resonance  of  a particular  transition,  the 
other  terms  make  less  than  1%  contribution  to  the  shift  and  we  may  accurately 
replace  the  resonant  contribution  to  6v  with  the  result  from  the  2x2 
degenerate  matrix  calculation  (Bonch-Bruevich,  1967), 


where  we  define 


(6) 


6 ••  O'  - A 


(7) 


5 


n'  = A'^  + Q.  ^ 

nm 

and 


(8) 


A = 


- V 


(9) 


The  sign  of  the  square  root  Is  taken  as  the  sign  of  A.  is  the  of f=resonance 

Rabi  frequency  of  population  inversion.  Note  that  the  maximum  shift  occurs 

at  A = 0 and  is  + /2.  Thus  very  close  to  the  n-m*"^  resonance  we  have 

— nm 


6v 

n 


1 

+ i ^ m n n m 


nm  2 I / / X 2 2 

m (v  -V  ,)  - V 

n m 


(10) 


with  the  dominant  term  6 having  the  resonance  limit  +(n  )/2  and  the 

remaining  terms  contributing  a non-resonant  background.  The  full  width 
in  A for  reaching  half  the  maximum  resonance  shift  is  easily  shown  to  be 
3n/2.  Thus  the  resonance  region  covers  the  range  in  A from  -3/4n  to  +3n/4. 

Figure  1 illustrates  on  an  energy  level  diagram  the  approach  to  resonance 
when  the  shift  in  level  n is  dominated  by  resonance  with  level  m as  A 
goes  to  zero.  6 and  are  as  defined  in  Eqs.  (7)  and  (8).  As  drawn, 

V is  less  than  the  transition  frequency  that  A is  positive. 

The  nth  level  is  repelled  and  pushed  down  by  5/2  while  the  mth  level  is 
repelled  and  pushed  up  by  6/2.  Photon  levels  at  + A and 

also  displaced  by  6/2.  Thus  the  pair  of  terms  of  the  irradiated  molecule  has  been 

split  into  four  terms  by  the  presence  of  the  light  field  and  there  are  four 

allowed  transitions  possible  as  shown:  two  at  the  Irradiation  frequency  v 

and  two  at  v + ft*  . As  A goes  to  zero,  each  term  is  split  by  6 •* 

and  the  transition  frequencies  are  v,  v + ft  . In  the  resonance  region 

— nm 

the  split  terms  represent  "dressed  atoms",  l.e.  a mixture  of  the  molecular 
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state  and  the  radiation  field.  The  molecular  part  of  the  split  term  changes 

from  being  predominately  in  the  down  (up)  shifted  part  of  the  lower  (upper)  term  as 

V approaches  from  below.  At  resonance  the  molecular  part  of  the 

state  is  equally  in  the  split  terms.  As  v goes  above  molecular 

part  of  the  state  shifts  rapidly  to  predominately  in  the  up  (down)  shifted 

part  of  the  lower  (upper)  term.  The  width  in  A of  this  crossover  is 

+(.31)f2  around  A=0.  This  is  discussed  in  Courtens  and  Szoke  (1977). 

— mn 

The  above  light  induced  terms  can  be  seen  in  spontaneous  emission  as 
a triplet  (Grove,  1977)  and  at  low  fields  with  a tuned  probe  beam  in  absorption 
as  a triplet  (Wu,  1977).  The  term  splittings  could  be  individually  studied 
by  observing  in  emission  or  absorption  a transition  which  began  or  ended  on 
one  of  the  resonantly  split  terms,  in  which  case  a doublet  would  be  observed. 

The  term  splitting  observation  was  anticipated  in  a calculation  by  Mollow 
(1972)  and  a first  laser  experiment  reported  in  the  photoionization  of  Na 
by  Moody  and  Lambropoulos  (1977). 

For  convenience  we  have  collected  in  Appendix  A some  numerical  conversion 

factors  useful  for  laser  induced  Stark  shift  calculations.  In 

this  report  the  key  to  the  frequency  scale  of  shifts  is  to  note  that  a 
2 8 2 

100  MW/cm  = 10  Watt/cm  incident  intensity  beam  corresponds  to  an  electric 
field  amplitude  of  275  KV/cra  or  915  Stat  Volt/cm.  The  Rabi  frequency  for  such 
a field  for  a 1 Debye  dipole  matrix  element  is  A. 60  cm 

III.  Application  of  the  Theory  of  Stark  Shifts  to  Diatomic  Molecules. 

In  order  to  apply  the  theory  to  any  species  we  must  reduce  the  term 
(energy  levels)  values,  the  selection  rules,  and  the  transition  intensity 
data  to  dipole  matrix  elements.  The  data  is  not  uniformly  available  for 
molecules  but  for  the  cases  of  the  diatomic  molecules  used  in  lasers  it  is 
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possible  to  collect  the  necessary  constants  from  the  literature  and  so  to 
evaluate  the  Stark  shifts  of  Eqs.  (5)  and  (8)  for  Illustrative  specific  cases. 

It  is  traditional  to  follow  the  notation  of  Herzberg  (1950)  for  the 
molecular  spectral  constants.  We  here  consider  diatomic  molecules  in  their 
ground  electronic  states  for  which  the  vibration-rotational  terms  are  labeled  by 
vibrational  quantum  number  v and  rotational  (angular  momentum)  quantum 
number  J.  The  standard  expressions  for  the  v , J labeled  terms  are  given 
in  Appendix  B together  with  the  spectroscopic  constants  used  in  this  report. 

We  note  here  that  these  terms  are  degenerate  in  magnetic  quantum  number 
Mj  = M.  As  we  discuss  below  this  degeneracy  is  removed  in  applied  electric 
fields. 

We  assume  that  the  selection  rules  for  diatomic  molecule  transitions 
are  the  same  as  the  selection  rules  for  dipole  allowed  transitions  of  a 
rotating  harmonic  oscillator  (Herzberg,  1950) : 


Av  = 0,  + 1 
AJ  = + 1 
AM  = 0,  + 1 


Thus  each  v,  J term  is  coupled  to  six  other  terms  by  dipole  allowed 
transitions. 

The  Av  = 0 transitions  correspond  to  pure  rotational  transitions. 
Though  they  have  large  transition  moments  from  the  static  dipole  moment  of 
the  molecule,  they  are  far  off  resonance  for  optical  radiation  at  frequencies 
near  the  vibration  rotation  frequencies.  Using  Eq.  (5)  with  v >>  I rot: . 

we  have 


(11) 


ll, 


I 


f 

h 


with 


and  p= 
aider  DF 
this  case 


■ i <Ef>' 


the  static  dipole  moment  of  the  molecule. 


with  = 1.82  D,  in  light  of  intensity  of 
n = 4.83  cm  ^ and  6v  = 2.8  x 10  ^ cm  ^ for 


As  an  examp Ir  r on 
10®  W/CB^.  Ic 

V • Ul  “ r* 


e 

This  shift  is  independent  of  J and  weakly  dependent  on  v.  It  is  muct. 
smaller  than  the  Doppler  width  of  the  vib.-rot.  transitions  in  DF.  TYv* 
molecules  HF  and  CO  give  smaller  shifts  and  so  the  pure  rotational  (.'.v  • 0 
contribution  to  the  Stark  shifts  will  be  neglected  in  the  following  taitu- 
latlons. 

The  remaining  four  transitions  which  couple  terms  to  the  v,  J ter* 
are  illustrated  in  Fig.  2 and  numbered  1 to  4.  In  the  conventional  notar  ■ r 
the  transitions  are  labeled  P if  the  upper  state  has  J decreased  b>  ; 
and  R if  the  upper  state  has  J increased  by  1 and  the  lower  state  ; 
is  given.  We  augment  this  notation  by  adding  the  lower  state  v to  i-**:  ■ 
the  band.  This  gives  the  notations: 


i 

v’{^}  (J") 

1 

“ VP(J) 

2 

vR(J) 

3 

(v-l)P(J+l) 

4 

(v-l)R(J-l) 

We  can  label  the  transitions  in  our  calculations  with  1*1  to  * smd 
term  differences  Av^  where 

Av,  - V -v^  " T(v",J")  - T(v  ,J  ). 

1 n m 11 


(12) 
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and  Rabl  frequencies  with  the  label  n referring  to  the  v,  J 

term  and  m the  appropriate  term  for  i = 1 to  4.  The  magnitudes  of  the 

transition  frequencies  order  as  < Av^  < IAV2I  < Av^.  Thus  the  general 

case  for  diatomic  molecules  in  the  vlb.-rot.  band  region  takes  the  form  of 

Eq.  (5)  with  just  four  terms  and  a fixed  order  of  transitions. 

An  example  of  the  full  range  of  term  shifts  is  shown  for  the  term 

v"  = 2,  J"  = 10  of  DF  in  Fig.  3 as  a function  of  the  applied  (laser) 

3 2 

radiation  frequency  at  an  intensity  of  10  W/cm  . Notice  that  this  plot  combines 
the  sign  of  the  shift  with  a log  scale  by  overlapping  the  6v  scales  at  ‘ 

+ 10  ^ cm  Fig.  3 shows  that  the  resonant  shifts  are  > .1  cm  ^ going  from 
small  negative  values  for  low  v to  small  positive  values  for  large  v > with 
sign  changes  at  each  resonance.  The  pattern  shown  in  Fig.  3 will  be  the 
same  for  all  diatomic  molecules:  only  the  specific  frequencies  and  on-resonance 
shift  magnitudes  will  change  with  molecule  and  term  considered.  We  will  return 
to  the  discussion  of  the  meaning  of  the  sign  changes  in  the  resonance  region 
after  we  take  up  the  question  of  the  proper  choice  and  calculation  of  transi- 
tion matrix  elements. 

IV.  Calculation  of  the  Transition  Dipole  Matrix  Element 

To  calculate  the  Stark  shifts  we  need  to  know  the  individual 
transition  dipole  matrix  elements  for  the  four  vlb.-rot.  transitions  which 
couple  to  each  term  v,  J.  Specifically,  we  need  the  squared  matrix  elements 
denoted : 

< v,J,MlM(r)-e^lv+l,  J+1,  (W^O.+l)  > . 

In  the  calculations  here  we  take  the  specific  case  of  linearly  polarized 

light  with  e parallel  to  z . This  corresponds  to  the  matrix  elements 

o i 


I 

A 
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< v.J.Mju^lv+l,  J+1,  M >lc^l  , 

since  the  z component  of  the  dipole  operator  gives  non-zero  matrix 
elements  only  for  AM  = 0.  The  angular  itegrals  of  the  rotational  states  are 
known  and  can  be  done  explicitly  (Townes  (1955)  and  Fenner  (1959))  leaving 
only  a radial  integral  as  a common  factor.  Thus  one  writes 

|<v,J,m|ij^1v+1,  J+1,  M >1^  = G(m,M)|<  v,j|u(r) |v+l,  J'  >|^  (13) 


where 


and 


(14) 


(15) 


The  matrix  elements  < v,j|p(r)|v',  J'  > are  called  the  radial  dipole  matrix 
elements  and  come  from  electronic  wavefunctlons  and  the  average  dipole  moment 
function  (Meredith,  1973).  The  radial  dipole  matrix  elements  are  further 
reduced  to  rotationless  dipole  matrix  elements  < v|p(r)|v'  > which  are 
Independent  of  J,  J' . It  is  the  rotationless  dipole  matrix  elements  which 
are  usually  tabulated  in  the  literature.  This  is  done  by  writing 

|<  v,Jlp(r) |v'J' ' >1 ^ = |<  v|p(r)|v'  >|^F  , (16) 


where  F if  a function  of  J which  is  nearly  unity  for  small  J.  In  the 
present  work  the  Oppenheimer  (1926)  result  is  used: 


F 


1 - Aym  [I-'Iyto-'Iy] 


(17) 
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n 

il 


where  y = (2Be/ue)  and  m given  in  Eq.  (15).  Discussion  of  improved 
F functions  can  be  found  in  Meredith  (1972)  and  (1973).  These  corrections 
are  less  Important  for  DF  and  CO  and  so  for  uniformity  they  have  not  been 
used  here.  Note  that  the  Oppenheiraer  F function  gives  increased  P branch 
matrix  elements  and  decreased  R branch  matrix  elements  as  J increases. 

For  high  J transitions  which  occur  in  diatomic  lasers  the  F factor 
gives  significant  corrections. 

Our  experimental  knowledge  of  dipole  transition  matrix  elements  usually 
comes  from  absorption  measurements  (Penner,  1959).  These  are  related  to 
Integrals  across  the  line  profiles  called  line  strengths  and  are  propor- 
tional through  known  factors  to  the  squared  radial  dipole  matrix  elements. 

Intensity  measurements  are  often  converted  to  spontaneous  emission  values 
called  Einstein  A coefficients  through  the  use  of  the  expression  (Penner, 
1959)  for  v'J'  ■>  vJ 


A(v'J' ,vJ) 


64Tr^v^|ml 

3h(2J'+l) 


1<  v,Jly|v',J’>| 


with  units  of  (molecule-sec.)  . Though  experimental  absorptions  may  be 


reduced  to  A coefficients  like  A(vJ,v'J')  the  theoretical  comparisons 


and  tables  of  coefficients  are  in  terms  of  rotationless  Einstein  coefficients 
A(v',v)  which  correspond  to  Eq.  (18)  for  the  case  |m|  = 1 (P(l)  Branch, 


J •=  1) . The  rotationless  matrix  elements  do  not  correspond  to  any  real 
transition  but  to  the  extrapolation  of  <v,j|vi|v' ,J'>  versus  m to  m = 0 
or  taking  vibrationless  matrix  elements  to  be  equal  to  <vj|vi|v'J'>  at  low 


J (J  = 0 or  1,  R(0)  or  P(l)). 

In  the  present  paper  the  needed  dipole  matrix  elements  were  calculated 
from  tables  of  rotationless  Einstein  coefficients  for  HF,  DF  and  CO 
found  in  Meredith  (1973)  and  Herbelln  (1974a,  1974b).  The  Einstein  coefficients 
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used  are  collected  in  Table  BII.  of  Appendix  B.  The  proceedure  used  was: 

(a)  Calculate  the  rotationless  dipole  matrix  elements  from  the  rotation- 
less A coefficients  using 

<vlu|v'>  = ^ A(v’,v)  . (19) 

641:  V 

(b)  Calculate  the  full  dipole  matrix  element  from  (a)  using  Eq.  (16) 
with  the  appropriate  F factor  from  Eq.  (17). 

(c)  Calculate  the  z component  dipole  matrix  elements  from  (h)  using 
Eq.  (13)  with  the  appropriate  G(m,M)  from  Eq.  (14). 

As  a check  on  the  results  of  steps  (a)  and  (b)  of  this  proceedure,  I 
calculated  the  A(vJ,v’,J’)  for  the  OP(IO)  transition  in  DF  and  found 
33.7  sec.  ^ compared  to  the  value  measured  by  Bonczyk  (1975)  of  34  ± 4 sec. 
Also  I calculated  <vj|pjv'J’>  for  the  4P(10)  transition  of  CO  and  found 
a value  of  .239  D compared  to  .238  ± .003  D measured  by  Weisbach  and 
Chackerlan  (1973). 

For  a given  (v'-v")  band  and  line  P(J")  or  R(J")  the  relative 
intensity  of  the  Stark  split  M M components  will  be  determined  by  the 
function  G(m,M)  of  Eq.  (14).  The  contributions  to  the  Stark  shift  of  the 
v",J"  term  will  also  be  proportional  to  G(m,M)  when  v is  off-resonance 
so  that  6v  is  proportional  to  the  squared  matrix  elements.  When  v is 
on-resonance  the  shift  of  each  M labeled  term  will  be  proportional  to 
G(m,M)**  . These  two  cases  are  illustrated  for  a P(10)  transition  in  Fig.  4 
(non-resonant)  and  Fig.  5 (resonant).  The  lifting  of  M degeneracy  is  well 
known  in  microwave  Stark  measurements  (Townes,  1955)  and  was  included  in 
detail  in  the  lineshape  calculations  of  Skribanowitz  et  al.  (1972).  Note 
that  this  degeneracy  lifting  occurs  in  atoms  as  well  as  molecules  and  so 
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is  always  to  be  expected  in  the  Stark  lineshapes  unless  M state  selection 
has  occurred  before  the  species  are  illuminated  with  the  laser  field.  Its 
effect  on  laser  induced  Stark  shift  spectra  lineshapes  was  not  discussed 
except  in  the  work  of  Grove  et  al.  (1977). 

Since  the  first  detection  of  laser  induced  Stark  shifts  in  molecules 
will  probably  occur  in  a gas  inside  a laser  without  state  selection, 
the  calculations  here  have  been  done  with  the  maximum  AM  = 0 matrix  element 
corresponding  to  M = 0.  The  actual  lineshapes  will  tend  to  be 
triangular  as  shown  in  Figs.  4 and  5.  The  first  moments  of  these 

"lineshapes"  are  slowly  varying  with  J and  fall  at  0.82  of  the  maximum 
shift  for  the  non-resonant  case  and  0.90  of  the  maximum  shift  for  the  resonant 
case  for  J=5  to  15.  Thus,  for  the  high  J cases  important  in  diatomic  mole- 
cules at  resonance,  half  of  the  intensity  of  a Stark  shift  will  occur  between 
0.90  and  1.00  of  the  maximum  shifted  component  (M=0)  and  we  need  only  calculate 
the  M=0  component  to  estimate  the  spectrum  of  shifts. 

V.  Stark  Shift  Calculation  Results 

In  this  section  the  numerical  results  of  the  calculated  Stark  shifts 
are  presented.  Specifically,  Eqs.  (5)  and  (8)  and  used  with  M = AM  = 0 
dipole  matrix  elements  calculated  from  Einstein  A coefficients  as  described 
in  Section  IV.  The  overall  pattern  of  the  shifts  for  a term  (v,J)  has 
been  shown  in  Fig.  3 for  the  term  v = 2,  J = 10.  This  shows  the  four 
resonances  which  will  occur  at  the  four  transitions  which  couple  to  this 
term. 

Since  the  2P(10)  line  of  DF  is  a strong  laser  line,  we  have  the  near- 
resonance region  of  this  line  displayed  in  greater  detail  in  Fig.  6.  The 

line  occurs  at  2494.37  cm  The  shifts  in  the  (2,10)  term  displayed  on 

-3  -1 

the  vertical  axis  cover  the  range  ±(1  to  10  ) cm  . Notice  that  for  the 
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two  laser  intensities  assumed,  the  off-resonance  shifts  vary  as  the  intensity 

(lOx)  while  the  resonant  shift  vary  as  the  square  root  of  the  Intensity.  Also 

indicated  on  the  figure  is  the  Doppler  width  5v^  (full  width  at  half  maximum) 

for  Av  = 1 transitions  of  DF.  (Calculations  and  table  of  Doppler  widths 

8 2 

is  given  in  Appendix  C.)  One  can  see  that  for  10  W/cm  intensity  the 

Stark  shift  should  be  resolvable  for  v within  ±5  cm  ^ of  the  2P(10)  line.  At 
9 2 

10  W/cm  even  the  non-resonant  shift  is  greater  than  the  Doppler  width 

over  a range  of  ±20  cm  Figure  7 shows  the  (2,10)  term  shift  near  2P(10) 

8 2 

on  a linear  scale  In  the  range  where  the  shift  at  10  W/cm  is  greater 

than  the  Doppler  width.  The  on-resonance  term  shifts  are  equal  to  half  the 

"1  8 2 

respective  Rabi  frequencies  at  the  two  intensities:  .188  cm  at  10  W/cm 
and  .594  cm  ^ at  10^  W/cm^. 

For  other  terms  of  DF  and  other  molecules  the  shifts  are  similar. 

Table  I collects  the  important  parameters  for  the  lower  state  terms  of 
several  strong  laser  lines  in  the  molecules  HF,  DF  and  CO.  The  lines 
were  selected  from  the  description  of  laser  line  strengths  in  the  literature: 

DF  lines  from  Kwok  (1970)  and  CO  lines  from  Yardley  (1971).  The  lasing 
transition  corresponds  to  1=1  in  the  table  The  labels  I = 1 to  4 cor- 
responds to  the  numbering  of  the  transitions  in  Fig.  2 and  Section  III. 

The  matrix  elements  are  for  the  AM  = M = 0 transitions  corresponding  to 
the  maximum  shift  and  maximum  Intensity  transitions.  The  molecular  parameters 
used  in  these  calculations  are  given  in  Appendix  B.  From  this  table  one  sees 
that  the  DF  results  presented  graphically  in  Figs.  3,  6 and  7 are  typical. 

The  CO  molecule  has  larger  matrix  elements  and  Rabi  frequencies  but  only 
by  2x.  The  term  shift  at  resonance  with  Av^^  is  given  with  = ±Dj^/2 
and  the  non-resonance  shift  coming  from  the  three  term  sum  over  non-resonant 


contributions . 
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j.  Table  I . Calculated  resonant  case  Stark  shifts  for  particular  lines 

and  terms  of  HF,  DF  and  CO  assuming  a laser  intensity  of 
8 2 

10  watts/cm  = e = 275  KV/cm.  The  dipole  matrix  element 
y(I)  corresponds  to  = 0 ] | v' , J',n=0>  | with 

J"  = J 2’  Av  and  n are  given  in  cm  y in  esu. 
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Table  I (cont.) 
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VI.  Conclusions  ' 

When  this  work  was  proposed  we  were  concerned  with  the  possible  effects  the  ' 

then  estimated  large  Stark  shifts  might  have  on  the  laser  output  spectrum, 
laser  turn-on  transients  and  atmospheric  propagation  of  high  intensity  lasers 
such  as  some  of  the  chemical  lasers  which  run  on  diatomic  molecular  vibrational- 
rotational  transitions.  After  careful  study  of  the  theory  of  these  Stark 
shifts  (Hollow,  1972a,  1972b)  and  the  atomic  Na  experimental  verifications 
it  is  clear  that  even  though  the  terms  of  the  molecules  inside  a high  energy 

laser  are  surely  split,  the  lasing  emission  continues  to  occur  at  the  dominant  , 

component  of  the  emission  of  the  coupled  terms  which  is  at  the  natural  fre- 
quency of  the  transition.  There  is  spontaneous  emission  which  will  be  occuring  j 

as  side  light  with  the  Stark  splittings  but  the  side  bands  are  not  able  to 
compete  in  stimulated  emission  with  the  stronger,  unshifted  component.  On 
the  propagation  question,  careful  study  of  the  sign  of  the  shift  given  by 
the  general  equation,  Eq.  (5),  shows  that  the  terms  always  move  so  as  to  shift 
the  molecular  transition  away  from  the  incident  radiation  frequency.  Thus 
any  absorption  line  causing  beam  attenuation  in  atmospheric  propagation  will 
be  Stark  shifted  away  from  the  high  intensity  laser  frequency  leading  to  a 
drop  in  attenuation  as  long  as  some  form  of  heating  or  higher  order  process 
doesn't  change  the  population  in  the  initial  state  of  the  near  resonant  transition. 

The  area  of  current  technical  interest  which  may  find  fruitful  application 
of  the  present  calculations  is  that  of  multiphoton  ionization;  it  is  clear 
that  the  high  intensities  used  in  these  studies  are  capable  of  causing  large 
shifts  in  energy  levels  (terms)  of  the  molecule.  If  the  multistep  process 
depends  on  some  intermediate  step  resonances  then  the  probability  for  the 
excitation  or  ionization  can  be  critically  dependent  on  the  intensity  of  the 
beam  due  to  Stark  shifts.  Moody  and  Lambropoulos  (1977)  showed  a simple 
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example  for  multiphoton  Ionization  of  Na.  Earlier,  Liao  and  Bjorkholm  (1975) 
had  demonstrated,  at  low  power  and  high  resolution,  the  sharp  resonance  in 
the  intermediate  state  which  could  be  Stark  intensity  tuned  in  two  photon 
absorption.  Papers  have  begun  to  appear  on  coupled  one-  and  two-photon 
resonant  dynamic  Stark  splitting  (Wong,  1977). 

It  seems  that  a serious  present  limitation  on  our  understanding  occurs 
because  as  yet  there  is  no  demonstration  of  dynamic  Stark  shifts  of  laser 
irradiated  molecules.  It  is  hoped  that  the  present  calculations  might  serve 
to  encourage  some  careful  experiments  on  diatomic  molecules  where  the  spectral  / 

constants  are  well  enough  known  to  allow  quantitative  comparison  of  experiment 
and  theory.  It  seems  that  any  scheme  which  can  selectively  and  efficiently 
produce  desired  photoionization  will  need  to  consider  the  dynamic  Stark 
shifts  of  molecular  terms. 
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APPENDIX  A 

Numerical  Conversion  Factors 

A-I.  Intensity-field  Conversion 

For  a field  e = cos (kx-wt) , the  average  Poyntlng  vector  magnitude 
gives  the  average  Intensity  as 

S * |<^>|  = (ergs/cm^)  , (A.l) 

in  cgs  units.  Noting  that  1 Stat  Volu/cm  = 300  Volt/cm  we  find  the  useful 
relations: 

e^(v/cm)  = 27.46  i^(W/cm^)  , (A. 2) 

and 

e^(Stat  V/cm)  = .0915  /s(W/cm^)  . (A. 3) 

A-II.  Dipole  Units 

Traditionally  dipole  moments  and  dipole  matrix  elements  are  given 

”18 

In  units  of  Debye  (D) , where  1 D = 10  esu-cm. 

A-III.  Rabl  Frequencies 

In  the  cgs  system  the  product  of  v in  esu-cm  and  e In  Stat  V/cm 

Is  an  energy  In  ergs.  For  a field  of  915  Stat  Volt/cm,  corresponding  to  a 

8 2 “18 

light  beam  Intensity  of  10  W/cm  and  p of  1 D = 10  esu-cm,  pe  = 

“18  “1 
915  X 10  erg.  Expressed  as  a Rabl  frequency  in  cm  we  have  0^  = 

(pe)/(hc)  = 4.60  cm  This  result  is  useful  In  scaling  other  cases  with 

different  p in  Debye,  and  intensity  S since 


APPENDIX  B 


Spectroscopic  Notation  and  Constants  for  Diatomic  Molecules 


Following  Herzberg  (1950)  we  write  the  term  value  of  the  v - 

tional  quantum  number),  J (angular  momentum  quantum  number)  as 

T = G(v)  + F^(J)  , (h.. 


with 


G(v)  = w (v+^s)  - to  X (v44j) 
e e e 


2 


(h.: 


F^(J)  = B^J(J+1)  - D^J^(J+1)^  + 


( B.  • 


where 


B = B - a (vHJ-5)  + ... 
V 6 e 


(B.*. 


and 


with 


D = D + B (v-P-s)  + 

V e e 


D 

e 


4 B ^/(o  ^ 
e e 


(B.  5 I 


(B.6i 


(0  , (OX,  B , a are  tabulated  in  Hertzberg  (1950)  and  more  recently  In 
e 66  6 6 

the  A.I.P.  Handbook  (Gray,  1972).  3^  is  much  less  than  so  we  neglect 

it  and  all  the  higher  terms  in  anharmonlclty  and  vibration-rotation  coupling 
The  spectral  constants  used  in  this  paper  are  collected  in  Table  B-1  for  the 
molecules  HF,  DF  and  CO.  The  Einstein  coefficients  (rotationless)  for 
the  Av  = 1 transitions  used  in  this  paper  are  collection  in  Table  B-1 I. 


i A-4 

r 

i! 

il 

ll 

i: 

Table  BI.  Spectroscopic  parameters  used  in  calculations  of  terms''  \ 
all  in  cm  ^ units. 


HF 

DF 

CO 

u> 

e 

4138.32 

2998.19 

2169.8233 

U X 

e e 

89.88 

45.76 

13.2939 

B 

e 

20.9557 

11.000 

1.931271 

a 

e 

.798 

.2907 

.017513 

D 

e 

2.149(-3) 

5.923(-4) 

6.120(-6) 

(a)  From  Gray  (1972). 

(b)  . 


I 


A- 5 


Table  BII,  Rotationless  Einstein  Coefficients  A(v+l,v)  in  sec'^. 


V 

MOLECULE 

HF 

DF 

CO 

0 

189. 

54.4 

33.4 

1 

322.5 

97.3 

64.2 

2 

406.1 

129. 

92.5 

3 

446.3 

151. 

118. 

4 

449.8 

163. 

141. 

5 

423.6 

167. 

162. 

6 

374.5 

163. 

180. 

7 

310.0 

153. 

196. 

8 

237.4 

138. 

209. 

(a)  From  Meredith  (1973) . 

(b)  From  Herbelin  (1974  a) . 

(c)  From  Herbelin  (1974  b) . 
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Doppler  Linewidths 


The  Doppler  width,  full  width  at  half  maximum  (FWHM) , of  the  Gaussian 
lineshape  of  a gaseous  molecule  is  given  by 


Av 

V 

o 


2R  logg(2) 

2 

c 


(C.l) 


with  the  transition  frequency  , c the  velocity  of  light,  R the  molar 
gas  constant,  T the  absolute  temperature,  and  M the  molecular  weight. 

For  T in  Kelvin  and  M in  gm/mole,  we  have 


^ = 7.162  X 10~^  (C.2) 

o 

Table  C-I  gives  the  values  of  Doppler  widths  for  HF,  DF  and  CO  at 
= 0)^  corresponding  to  Av  = 1 transitions. 


j. 


Table  Cl.  Doppler  widths  at  T = 300  K,  Av  = 1 transitions. 


Molecule 

M (gm/mole) 

Vq  (cm  ) 

Av  (cm 

HF 

20 

4138. 

.011 

DF 

21 

2998. 

.008 

CO 

28 

2170. 

.005 

(a)  Av  = Vg (7. 162x10" FIVHM. 
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g.  2 Schematic  diagram  of  the  vibration-rotation  terms  coupled 

by  allowed  transitions  to  the  term  v,J.  The  numbering  inside 
circles  corresponds  to  the  running  index  i 1 to  4 in  the 


Fig.  5 Relative  intensity  vs.  relative  Stark  shifts  for  the  M corapon 
ents  of  a P(10)  line  AM=0  transitions  at  resonance.  Numbers 


solid  curve  for  10  W/cm  and  broken  curve  for  10  W/cm' 
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